#! /usr/bin/env python

import pyx
import pyfeyn
from pyfeyn.user import *

fd = FeynDiagram()
processOptions()

in1 = Point(-6,  2)
in2 = Point(-6, -2)
in11 = Point(-2,  0)

out11= Point(2,0)

out1 = Point(6, 2)
out2 = Point(6, -2)
out = Point(6, 0.5)

vtx = Vertex(3, 0.5, mark=CIRCLE)

Fermion(in1,in11).addArrow().addLabel(r"\Pq")
Fermion(in11,in2).addArrow().addLabel(r"\Paq")
Gluon(in11,out11).addLabel(r"\Pg")
Fermion(out11,out1).addArrow().addLabel(r"\Pqt/\Pqb")
Fermion(out2,out11).addArrow().addLabel(r"\Paqt/\Paqb")
Higgs(vtx,out).addLabel(r"\PHiggs")


fd.draw("./basoc_Higgs_production_Feynman2.pdf")

